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In this paper we apply the emerging- consensus understanding of the fermionic self energy deduced 
from angle resolved photoemisssion spectroscopy (ARPES) experiments to deduce the implications 
for orbital diamagnetism in the underdoped cuprates. Many theories using many different starting 
points have arrived at a broadened BCS-like form for the normal state self energy associated with 
a d-wave excitation gap, as is compatible with ARPES data. Establishing compatibility with the 
f-sum rules, we show how this self energy, along with the constraint that there is no Meissner effect 
in the normal phase are sufficient to deduce the orbital susceptibility. We conclude, moreover, that 
diamagnetism is large for a d-wave pseudogap. Our results should apply rather widely to many 
theories of the pseudogap, independent of the microscopic details. 

PACS numbers: 



INTRODUCTION 

The origin of the pseudogap in high temperature su- 
perconductors is still very much under debate. While 
there are indications of an alternative order parameter 
[J, Q it is widely accepted that the diamagnetism which 
is associated with the pseudogap is indicative of some 
form of precursor pairing. Nevertheless, even within the 
various schools which posit preformed pairs there is con- 
troversy about their nature and origin. What would be 
particularly useful, given this multitude of scenarios is to 
build an understanding of the pseudogap phase based on 
aspects of phenomenological (and microscopic) theories 
about which there is some degree of convergence. 

In this paper we focus on the diamagnetism of the 
cuprates (in the linear response regime) and show how 
it is connected to an emerging- consensus understanding 
of angle resolved photoemission experiments. These 
experiments have demonstrated that a simple broadened 
BCS form for the self energy in the presence of a pseu- 
dogap Apg_k 
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works rather well. Indeed, this form of the self energy 
is widely subscribed to in diverse theories of the 

pseudogap. Here we exploit this equation (and the re- 
lated fermiology constraining the bandstructure a) in 
conjunction with the important constraint that there is 
no Meissner effect in the pseudogap phase, to arrive at 
a form for the orbital susceptibility. This relation be- 
tween transport and self energy is well known from the 
Ward identities. In a closely related fashion we justify 
our transport expressions for the diamagnetic response 
by showing that the related current-current correlation 
function analytically satisfies the transverse and longitu- 



dinal f-sum rules. In this sense this understanding of dia- 
magnetism should have a generality which goes beyond 
a particular microscopic or phenomenological approach. 

Recent experiments measuring the orbital susceptibil- 
ity in a variety of cuprate superconductors report anoma- 
lously large diamagnetism [2, ;8j). This diamagnetism 
onset is well correlated with the onset of an enhanced 
Nernst signal above Tc in the same materials [9]. It is 
not likely that low dimensional, critical fluctuations are 
responsible, because the diamagnetism persists beyond 
the expected critical regime. Moreover, there are reports 
of strong non-linear effects, although they will not be the 
topic of this paper. 

This pseudogap phase has been very systematically 
studied in angle resolved photoemission spectroscopy 
(ARPES) experiments which lead to the phenomeno- 
logical expression in Eq. ([I} as well as the underlying 
fermiology. An important feature in the precursor super- 
conductivity approaches is the Fermi "arcs" , (as distin- 
guished from Fermi "pockets"). These are understood 
to arise from the phenomenological broadening factor 7 
in Eq. (jlj which leads to a smearing out of the d-wave 
nodes. Also relevant are the collapse of the arcs below 
Tc and the related "two gap" phenomena which appear 
in the superconducting phase. Importantly, they fol- 
low naturally [13, HI from Eq. ^ if one adds to Epg 
a co-existing conventional condensate contribution E^c, 
(of the same form but with the order parameter contri- 
bution, Agj, appearing and the associated 7 = 0). The- 
oretical approaches to diamagnetism have not incorpo- 
rated this ARPES constraint, largely because they are 
based on thermodynamic rather than transport schemes. 
Moreover, for the most part the emphasis has been on the 
contribution of bosonic pairs either in fluctuation [l^ - [l^ 
or real space pairing schemes. [l3|. Nevertheless, physi- 
cally, one might expect the total diamagnetic response to 
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reflect fermions as well as these pair correlations. Indeed, 
this has to be crucial when the fermionic self energy is 
incorporated and sum rules are addressed. Alternative 
theories based on vortex liquids [l6| or d-density wave 
states [17| have addressed complementary physics. 



Overview of Theory and Results 

It is useful to begin with an overview of the connection 
between ARPES and diamagnetism. The reasoning is 
relatively simple and is best communicated by presenting 
the theory first in summary form. 

For the sake of notational simplicity, the equations we 
present here are for the negligibly small 7 limit, although 
the general line of reasoning can readily accomodate any 
size 7, as discussed in the Appendix. The transport- 
based expression for the diamagnetic response is 



- lim Re 
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(2) 

Here x'^"' is written in terms of the current-current cor- 
relation function ^(q, w) and the diamagnetic current 
tensor ^ /m. The absence of a Meissner effect is equiv- 
alent to the vanishing of the electromagnetic response 
kernel above Tc at zero momentum and frequency, i.e. 
the response kernel must satisfy 



^(0,0) + — ^0,T >Tc 
m 



(3) 



When Eq|3]is satisfied, the first nonzero term in the nu- 
merator of Eql2] comes in at ©(q^), rendering x'*"* well- 
behaved. The analogue statement below Tc is that the 
right hand side of Eql3]equals ^^/m, the superfluid den- 
sity tensor. 

It is helpful to recast the expression for /m in a way 
that closely mirrors the structure of a response function: 

— ~ 2 > „, ^, Gk can be rewritten after some straight- 
en r)\rr)\r 
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forward algebra as 
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Now let us consider the relation to Landau diamagn- 
tism. Ignoring lattice effects for the moment, the stan- 
dard expression for the orbital susceptibility is of the form 
( 
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where n and m represent the density and mass of the 
charge carriers in Landau orbits. The source of large 



conventional fluctuation diamagnetism [12| is the char- 
acteristic size (r^) which is presumed equal to the cor- 
relation length and which diverges as the transition is 
approached. Consequently, even in the absence of fluc- 
tuations, one might expect that s-wave pairs would have 
significantly less diamagnetism than the more extended 
d-wave pairs of the cuprates, as we demonstrate here. 

We will show that the orbital susceptibility deriving 
from Eq. ([T]) arises from both fermionic and bosonic con- 
tributions. Diamagnetism comes predominantly from the 
latter and one can understand this physically as asso- 
ciated with the general enhancement of bosonic contri- 
butions to transport in the vicinity of Bose condensa- 
tion. A strongly peaked Bose distribution function (in 
momentum space) readily accomodates a redistribution 
of particles leading to large transport responses even in 
the presence of weak perturbing fields. Indeed, this is 
the origin of superconducting fluctuation contributions, 
in general. This is in contrast to fermionic transport, 
which is restricted by the Pauli principle. 

We note that [20| even in a zero gap normal state, 
bandstructure effects can yield a paramagnetic orbital 
susceptibility, particularly near Van Hove singularities. 
Thus one might expect the net diamagnetic contribution 
will be largest away from the Van Hove points. An ad- 
ditional effect of the fermionic contribution derives from 
the fact that an excitation gap will reduce the number 
of available fermions n. Because of nodal fermions, d- 
wave pairing in the pseudogap phase is expected to yield 
more diamagnetism. In this paper we will show how all of 
these effects combine to yield a rather large diamagnetic 
response strongly associated with a d-wave pseudogap. 

We present in this introduction an important inference 
which will be discussed in more detail later and which 
allows us to arrive at an extension of Eq|4] to finite mo- 
mentum and frequency in the form 
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One could anticipate such an answer using the standard 
normal state expression for the current-current correla- 
tion function (c.f. [^) or alternatively the counterpart 
for the standard BCS current-current correlation func- 
tion (where one has to be careful to enforce a Meissner 
effect, rather than its absence in the pseudogap phase). 
But the strongest support for Eq. ([6]) is the demonstra- 
tion that it analytically satisfies the f-sum rules. These 
are discussed in more detail in Appendix A. Importantly 
Eq. (|6]) provides the input one needs to arrive at trans- 
port properties (including the complex conductivity) that 
arise from Ward identity compatibility with the ARPES- 
derived self energy. 

It is useful to end this introduction with a summary 
figure addressing the implications of the pseudogap as 
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Figure 1: Summary Figure showing the efTect of the pseudogap for s and d-wave lattices via plots as functions of band filling 
or chemical potential, /j.. (a) corresponds to the assumed density of states, (b) to the normal state susceptibility while (c) and 
(d) illustrate the the difference between the pseudogap and normal state x"^^^ normalized by the absolute value of x"*"* in the 
normal state at the same /i. In (c) the pseudogap is taken to have an s-wave form fac tor. (d) The analogous plot but for a 
d-wave order parameter. 



it appears in Eq. ([T]) on the orbital susceptibility. The 
four panels in this figure correspond to (a) the assumed 
density of states as a function of energy and (b) the coun- 
terpart normal state orbital susceptibility as a function of 
band filling. Note that the simple Landau diamagnetism 
of jellium can give rise to paramagnetism in a tight bind- 
ing band, particularly near Van Hove singularities. In (c) 
we indicate the change in the orbital susceptibility asso- 
ciated with the pseudogap as a function of band filling 
for s-wave pairing. The counterpart figure for the d-wave 
case is shown in (d) . A comparison shows that the d-wave 
pseudogap is associated with substantially enhanced dia- 
magnetism. This figure will be discussed in more detail 
in the context of our numerical results. 



THEORY OF DIAMAGNETIC SUSCEPTIBILITY 

We begin by rewriting the expression for ^ /m: 



= -2 



^ dk dk 



(7) 



following an integration by parts. Next we rewrite the 
derivative of the Green's function using the identity 



dGx/dk — —Gl^dGj/ /dk, where the derivative of the 
inverse Green's function is straightforward to evaluate in 
terms of the ARPES-derived Yipg^K, 
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From Eq|3l the current-current correlation function at 
q, ifi = necessarily yields EqlH 

The structure of Eql4] is intriguing: the first term 
in parenthesis is the usual electromagnetic response of 
fermionic quasiparticles, while the second term arises 
from the presence of pseudogap correlations and has the 
appearance of being a correction to the bare electromag- 
netic vertex. 

We then posit Eql6] as the natural extension to arbi- 
trary four-vector Q, and this is supported by the trans- 
verse and longitudinal f-sum rules. For the latter we have 
to prove the longitudinal and transverse f-sum rules. Pre- 
viously [H, [13] we have focused on the second of these, 
which will be summarized here in the Appendix. We now 
address the longitudinal sum rule 
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We define = /(i^k) as the fermi function, == 
Ew±hfq/2, = ■Ck+q/2, and A^g = Apg^k±q/2- The nor- 
mal state fermions have dispersion Ey^ = + A^g(T). 

The longitudinal component of XJJ is defined as Xjj = 
q • Xjj • q = q • xjj • q/?^- We need several simple 
relations 
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along with 
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Changing variables p — s> p + ^ for the first and p — > p — f 
for the second term, yields the desired result 



Explicit Calculation of Diamagnetic Susceptibility 
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The electromagnetic response kernel has been con- 
structed so that there is no Meissner effect above Tc, 
only the coefficient in the series expansion of Eq|6] is 
needed to calculate x'^'^. The calculation is lengthy but 
straightforward, with the result that x'^'^ can be written 
as the sum of two terms, x'^''^ = Xo'^ + Sx'^^^ where 
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Here we have separated terms so that the first is the 
contribution independent of derivatives of the gap form 
factor, while the second term introduces these derivative 
contributions, which are notably absent in the s-wave 
case. In the above equations we have = i (l 
and m^^ = d'^^k/dk'^. 



The expressions Ea ll6IIT7l are quite complicated, but 
general inferences can still be made on the expected be- 
havior of the orbital susceptibility. We find that to a 
good approximation Eq |16l can be written as 
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and EqlT7]is well approximated as 

'l /90 
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It is seen from a numerical analysis that the largest con- 
tributing term to this last equation is the {d4>k/dky)'^ 
piece. 

NUMERICAL RESULTS 



The sum of Equations ((18)) and (|T9|) provides a rea- 
sonable approximation to the total orbital susceptibility 
Eq. (|16 p - (fT7|) . Moreover all that is needed is a the tem- 
perature dependent pseudogap parameter (and the chem- 
ical potential). It is our intention to first present results 



that are independent of microscopic details so the plots 
are representative of general theories in which we treat 
Apg and fj. as variables. We will simultaneously consider 
s and c?-wave pairing in both jellium and tight binding 
lattice models. For the latter we consider the simplest 
dispersion 
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-2t{cos{kx) + cos{ky)) — At' cos(fca;) cos(A;j,) — /(20) 



along with more realistic fits to the cuprate bandstruc- 
ture via ARPES data. 



Jellium models with s and d-wave pairing 

We begin by studying a jellium dispersion and an s- 
wave gap which avoids the complications of lattice effects 
and gap nodes. In Figl2l^a) we plot the orbital suscep- 
tibility from Eql^l Curves are shown for Apg = 0, 0.6, 
and 1.2Ep and normalized by xns, defined as the abso- 
lute value of the normal state orbital susceptibility when 
/i = 0. The system exhibits diamagnetism even in the 
normal state, reflecting the well-known Landau orbital 
susceptibility corresponding to minus one third of the 
Pauli susceptibility. The degree of diamagnetism is in- 
creased modestly as the gap size increases; this numeri- 
cally confirms the intuition that the bosonic term in EqlB] 
results in an enhancement of x^^'^- 

For the case of a c?-wave gap with a jellium dispersion 
we use Apg k = Apg cos 20 where (f) is the azimuthal an- 
gle. The resulting plot is shown in Figl^Jb). The curves 
are again displayed for Apg = 0, 0.6, and 1.2£^f and nor- 
malized by xns- In contrast to the s-wave case, where 
X grows as large as x'^''^ ^ "i-^XNS, x'^^'^ ^T-^Xns in 
the d-wave case, is roughly an order of magnitude larger 
than the normal state orbital susceptibility. 

We understand this contrast between the s and d- 
wave order parameters as physically resulting from the 
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Figure 2: The orbital susceptibility for a jellium dispersion 
with an s-wave (upper) and d-wave gap (lower). We consider 
Apg = 0,0.6, and 1.2Ef. The curves are normalized by Xns, 
which is defined as the magnitude of the normal state value 
of x'^'"' at M = 0. There is a finite orbital susceptibility even 
in the normal state. The diamagnetism is enhanced for all 
fi as size of the gap is increased. The d-wave gap symmetry 
results in a larger x'*'" than in the s-wave case. 
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Figure 3: Change in the orbital susceptibility with varying 
pseudogap for a jellium dispersion. Plotted is the relative 
difference between the pseudogap and normal state x"^'^ ^ 
a function of chemical potential /i for A = 0.6, and 1.2Ef- 
Upper panel is for an s-wave and lower panel for a d-wave 
gap. The presence of a pseudogap results in an increase of 
X'*'^ for all values of /i. The trends as a function of /i and A 
are the same but the overall magnitude of the enhancement 
is more significant for the d-wave case. 



increased pair size of the latter 2J, |25| , as reflected in 
the generic expression for orbital diamagnetism (Eql5|). 
The increased pair size reflects the nodal structure of the 
gaps; the correlation length of the pairs (which is not 
necessarily the same order of magnitude as the pair size) 
has a similar increase due to the nodes. It can be argued 
that the correlation length ^ oc A^^ where l^ps is the 
average of the gap along the Fermi surface. The quantity 
^FS will always be less relative to an s-wave order pa- 
rameter as long as the gap maximum is set to the same 
value as the s-wave gap that it is being compared against. 

Fig 13] shows the change in the orbital susceptibility as- 
sociated with a pseudogap for the case of jellium. Here we 
plot 5xNS = x(m) - XNs{^J') normalized by \xNs{^J')\ for 
the gap values A^^ = 0.6, and \.2Ep^ where the normal 



state susceptibility is at /i = 0. We consider A^ 



0.6 



and \.2Ep for the blue and black-dashed curves respec- 
tively. FiglSja) shows 5xns for an s-wave order param- 
eter and FiglHJb) plots the counterpart figure, 5xns, for 
the c?-wave case. For both the s and d wave cases when 
the filling is decreased (so the the normal state Landau 
orbital susceptibility is smallest) then the relative effect 
of the pseudogap appears more prominently. As antic- 
ipated, the overall diamagnetic contribution due to the 
d-wave order parameter is significantly larger than the 
contribution occuring for an s-wave pseudogap. 

In summary, for a jellium dispersion, we have seen 
that (i) the presence of a pseudogap leads to an enhance- 
ment of the normal state diamagnetism for most values 
of the chemical potential /i. Critically, (ii) a d-wave order 
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Figure 4: The orbital susceptibility integrand for a tight- 
binding dispersion at an angle 7ir/36 in the kx-ky plane, 
(top) The integrand for an s-wave order parameter of value 
Apg — 30, 60, and 90meV. The curves are normalized by A, 
defined such that the area under the Apg — 30 meV curve is 
unity. The curve for all values of Apg results in a paramag- 
netic contribution. The total area normalized by A of each 
curve is written alongside of it; the area takes the values 1, 
0.97, and 0.96 for Apg — 30, 60,90meV respectively, (bottom) 
The same plots for a d-wave order parameter. The curve for 
Apg — 60 and 90 meV results in an overall diamagnetic con- 
tribution. The total area normalized by A of each curve is 
written alongside of it; the area takes the values 1, —5, and 
— 10 for Apg = 30,60,90meV respectively. 



parameter results in a significantly larger diamagnetism 
than in the case of an s-wave order parameter. 



Tight-binding lattice 

It is useful to understand the integrand before address- 
ing the entire expression for the orbital susceptibility in 
a tight binding lattice. Figure U] plots the integrand asso- 
ciated with Ea ll6lfT71 as a function of radial distance in 
plane (normalized by the lattice spacing a) . We 
consider three different values of the gap (Apg ~ 30, 60, 
and 9QmeV) . Each of the curves is normalized by a num- 
ber A selected so that the curve corresponding to the 
smallest gap Apg = 30meV has area unity; the area of 
the curves is noted beside each. All curves correspond to 
a cut at 35 degrees above the axis. Fig|4lja) shows the 
integrand corresponding to an s-wave gap and FigU^b) 
corresponds to the d-wave case. 

For the former these cuts result in paramagnetic contri- 
butions to x'^"^ for all values of Apg and the extent of the 
contribution is largely unchanged as Apg increases: de- 
spite a tripling of the gap size, the displayed areas change 
by no more than 4% and lead to a paramagnetic contri- 
bution. 

The corresponding curves for a d-wave order param- 
eter are displayed in Fig. SJb). Although the form of 
the curve does not change qualitatively as the gap size is 
changed, the initial area of the curve is sufficiently small 
such that minor changes to the cuts result in significant 
changes to the nature of the overall contribution to x'^'^- 
as the gap increases from 30 to 90meV, the overall con- 
tribution changes from being paramagnetic to being 10 
times as large and diamagnetic. This trend is driven by 
the d-wave gap symmetry. 

We now return to FiglT] which plots the density of 
states versus energy (a) and the orbital susceptibilities 
in the gapless normal state (b) and in the presence (c) 
of an s- or(d) a d wave pairing. The last three panels 
are plots as a function of band filling or chemical poten- 
tial and the units arc in terms of the nearest neigh- 
bor coupling which is taken to be 300 meV. Panel (b) 
shows that when is large, the density of holes or elec- 
trons is sufficiently small such that the system behaves 
similarly to the jellium case and x'^^^ leads to the usual 
Landau orbital susceptibility result. There is a competi- 
tion [2^ , however, between the Ep and t energy scales, 
so that once lEpl/t < 1.8 the system becomes entirely 
paramagnetic. The term that contributes to this param- 
agnetism is weighted by the density of states evaluated at 
the Fermi surface and so it is dramatically enhanced by 
the Van-Hove point at Ep — 0. It is clear from the figure 
that the pseudogap-enhanced paramagnetism is substan- 
tial but only for the d-wave case. 

FiglSja) presents a plot of the orbital susceptibility for 
an s-wave gap and a tight-binding dispersion normalized 
by XNS, which is defined as in previous figures. The 
four curves correspond to the gap values Apg = 0, 30, 60, 
and 90ineV. The lattice dispersion Eql^D] is used with 
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Figure 5: (left) The orbital susceptibility for a tight-binding 
dispersion with an s-wave order parameter. The normal- 
ization xns is defined as the value of x"^'^ in the normal 
state when fi = 0. The curves correspond to gap sizes 
Apg = 0, 30, 60, 90 meV. There is very little change of x'^''' 
regardless of the gap size, (right) The orbital susceptibility 
for a tight-binding dispersion with a d-wave order parameter 
normalized by Xns- The d-wave form factor and the result- 
ing enhancement in the size of the pairs leads to a significant 
enhancement of the diamagnetism for ranges of fi away from 
the Van-Ifove point as the gap size increases. 



t = — 300mey and t' = 0. The s-wave pairing gap has 
little effect on x'^"^, with all four curves tending to overlap 
except near the Van-Hove point at /i w 0. The curves 
are entirely paramagnetic for a significant range of filling 
and the diamagnetic features at other filling are relatively 
small: an s-wave gap on a lattice is not enough to capture 
the observed experimental effects. 

The corresponding figure for a c?-wave gap is displayed 
in FiglUb). Gap and dispersion parameters that are 
identical to the s-wave parameters are used here. It 
is immediately seen that there is a significant enhance- 
ment of the diamgnetic susceptibility as the magnitude 
of the d-wave order parameter is increased. The Van- 
Hove point suppresses the extent of the diamagnetism 
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Figure 6: The behavor of the diamagnetism in a cuprate-like 
bandstructure. (upper panel) The model density of states 
showing that the Van-Hove point occurs at negative energies. 
Optimally doped samples have chemical potentials closer to 
the Van-Hove point than underdoped samples, as indicated 
by the arrows, (lower panel) The diamagnetic ssuceptibility 
as a function of /i for Apg = 0, 30, 60, and QOmeV normalized 
by the value of the normal state x'*'^ at /i = 0. 



and results in paramagnetism for small gap values and 
fillings near /i = 0; nevertheless, the diamagnetism per- 
sists over a wide range of filling for each of the gap val- 
ues considered here. For fillings that are near mid-band 
(i.e. fi « 500mel^), it is seen that there is a substantial 
increase in x*^'^ as compared to the normal state once 
a finite Apg is considered. The extent of this enhance- 
ment increases as one moves away from the Van-Hove 
point, mirroring the experimentally observed trend of the 
anomalous diamagnetism having the greatest impact in 
underdoped samples. 

We turn now to more realistic parameters for the band- 
structure of the cuprates in FigE] Here we use the 
ARPES-derived parameters ^ for Bi2212. With next 
nearest neighbor contributions one sees the broad trends 
are consistent with the earlier results. Note that the finite 
t' breaks particle hole symmetry and shifts the Van-Hove 
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Figure 7: The diamagnetism as a function of temperature 
for different "hole concentrations" parameterized via T* /Tc- 
The dashed Une shows the normal state value of t he sus- 
ceptibility. This diamagnetism is to be associated with the 
large transport response from pairs or bosons near conden- 
sation. We contrast the present precursor superconductivity 
scenario with the more more widely discussed phase fluctu- 
ation approach where the diamagnetism is associated with 
mesoscopic supercurrents. 



point in the density of states. 

FiglHJa) indicates a typical density of states for a 
cuprate system, along with the estimated positions for 
the chemical potentials. Here the Van-Hove point is lo- 
cated at negative energies. The general trend is that opti- 
mally doped systems will be closer to the Van- Hove point 
than underdoped systems ( [26| ) . Interestingly, this along 
with the stronger pseudogap in the underdoped case 
leads us to anticipate that underdoping enhances diamag- 
netism. The resulting orbital susceptibility is plotted in 
FigIMb) as a function of band filling for A = 0, 30, 60, and 
90meV. The curves are qualitatively similar to previous 
figures with the normal state x'*"' having both regions of 
paramagnetism and diamagnetism and becoming entirely 
diamagnetic for a sufficiently large value of Ap^. 



Temperature Dependence of Diamagnetism 

The above discussion has been quite general; we have 
essentially explored the consequences of Eq. ^ in all 
its generality. One can then inquire as to how this ap- 
plies to the high temperature superconductors. This re- 
quires that we establish the parameters Apg(T) and the 
fermionic chemical potential, n{T). Here, for definite- 
ness we use our preformed pair scenario [l^] which is 
based on stronger than BCS attractive interactions (con- 
sistent with small pair size and anomalously high pairing 
onset temperature T*). Once the pseudogap and chem- 
ical potential parameters are self consistently obtained 
27| . one accomodates a variety of dopings, by effectively 
fitting 27| the attractive interaction to match T* and 



Tc- For definiteness, we presume the band dispersion is 
associated with t = SOOmeT^ and /x = — 1.75t, chosen 
somewhat away from the Van Hove point. 

We present temperature dependent plots of the dia- 
magnetic response x'^'^'i i^i Figure [7] for four different 
dopings. Each curve is normalized by its normal state 
value, which is separately plotted as a dashed line. In- 
dependently of the particular parameters that are used, 
it is seen that the magnitude of x"^''' is enhanced even at 
temperatures well above T^. Importantly, this diamag- 
netism is not restricted to two dimensional models, as in 
fluctuation theories. Experiments as well re-enforce three 
dimensional critical behavior 28- 30|. 

In this way the present physical picture contrasts with 
the traditional fluctuation approach, in which one might 
expect large diamagnetism but only in the narrow criti- 
cal regime. Here, it is the stronger than BCS attraction 
which stabilizes these pair degrees of freedom (up to high 
temperatures T k T*) rather than the low dimensional- 
ity [12]. Importantly, in the present theory we compute 
the total conductivity based on the pseudogap self en- 
ergy, not the fluctuation corrections and in this way are 
analytically able to establish compatibility with the con- 
ductivity sum rules. 



COMPARISON BETWEEN DIAMAGNETISM 
AND CONDUCTIVITY 

The orbital susceptibility is not the only transport 
property which can be associated with the pseudogap self 
energy of Eq. ([TJ. We have previously discussed the op- 
tical and THz conductivity 22 , 23| . Here we concentrate 
on the challenges raised by recent experimental papers 
[31I [3^ which have pointed out the seemingly contra- 
dictory behavior implicit in the dissipative conductivity 
and the orbital susceptibility. The authors of Ref. |31 
deduce "However, if the diamagnetism signal above Tc 
is solely due to superconducting correlations then it is 
a well posed challenge to explain the lack of straightfor- 
ward correspondence to conductivity." In this section, we 
address this challenge by looking simultaneously at the 
linear diamagnetic response and the dissipative conduc- 
tivity. 

In this regard, it is important here to incorporate the 
ordered phase. We can anticipate that we now have two 
distinct contributions to the fermionic self energy S(k, uj) 



S(k,w) 



-17 - 



pg,k 



A2 



6 



«7 



6 



(21) 



The first of these is associated with the normal state and 
the second with the condensate. 

Following the analysis of this paper (above Tc) we may 
similarly arrive at the below Tc counterpart directly from 



10 




Figure 8: Numerically generated plots of (a) the imaginary conductivity 02 as a function of temperature at constant energy 
uj. The curves are normalized by ctq, the value of 02 at zero temperature with uj — ImeV. (b) The imaginary conductivity 
CT2 as a function of frequency for constant temperature T. The curves are again normalized by (Tq. (c) The conductivity and 
diamagnetic susceptibility as functions of Apg at ksT = 15meV. 



this self energy. Here for simplicity we rewrite P{Q) in 
the regime of very weak dissipation (7 w 0) and s-wave 
pairing where the behavior is more physically transpar- 
ent. 



k 



v-^ kk 



m 



E++E_ 



E+E. 



(1-/+-/-) 



E-\-E— 



E^ 



+ 



E-) 



E+E- 



E) 



f- 



(22) 



where /± - fiE±) and ^A^ = - A^^, e± = ?k±q/2, 
and E± — i?k±q/2 • The transverse sum rule can be shown 
to be precisely satisified, below as well as above Tc- The 
longitudinal sum rule cannot be easily proved below 
(due to collective mode effects) although it is analytically 
satisied in the normal state as shown earlier in this paper. 

Our interest here is on contrasting the diamagnetism 
with the dissipative conductivity 



(72 (w) = limRe 



(23) 



where the imaginary counterpart corresponds to the real 
part of the frequency dependent conductivity cti. Note 
that Eqs. ([2]) and ((23|) are superficially rather similar, but 
importantly very different. This difference derives from 
the different denominators and in a related fashion, the 
order of q, w limits. 

We stress that our below Tc extension can be derived 
3? 



microscopically [27|, 133 



However, the present more 
general approach based on the self energy indicates that 
the results should be rather generic. In support of this 
below- Tc generalization of the self energy, it is straight- 
forward to see from Eq. (|2T]) how the Fermi arcs will col- 
lapse to point nodes as soon as the system passes below 
Tc [l3| and that the general 2-gap phenomena, in which 
the nodal region is sensitive to A^c while the antinode 
region exhibits little temperature dependence fllj| are all 



direct consequences. In the same way one can explore the 
optical conductivity j22j and THz conductivity [23| , all of 
which are highly constrained by the simple self energy ex- 
pression in Eq. ([T|) and its extension below Tc in Eq. (|2ip . 
In addition, these expressions for the self energy lead to 
results for the specific heat and serve to constrain the 
tunneling characteristics, including quasi-particle inter- 
ference 136( 1 . Because the measured diamagnetism seems 
to be well correlated with the Nernst effect, it is impor- 
tant to note that a large Nernst response was addressed 
earlier and argued [37| to be associated with pre-formed 
pairs, as distinguished from normal state vortices. 

It is of interest then to compare the behavior of 
the conductivity and the diamagnetic response. FigH] 
displays our results for the cuprate models, based on 
Eqs. ^ and ([^5)) . FiglSJa) shows how the normal state 
cr2(w) behaves as a function of temperature for three dif- 
ferent low frequencies and FiglSj^b) plots the imaginary 
conductivity (T2 as a function of frequency for a range of 
different temperatures. At roughly Tc, we find that (T2 
shows a sharp upturn at low uj. Below Tc we find the 
expected supcrfluid frequency dependence <T2 cc Us/uj. 

We find that 172 (w) is modestly increasing with increas- 
ing w above Tc as seen in experiment, albeit slightly away 
from the transition. Experimental studies do reveal a 
small 10-15K range where a fluctuation contribution is 
visible. But except in this narrow temperature region, 
the observed behavior is not compatible with that ex- 
pected of a phase fluctuation contribution, where a uj~^ 
dependence would occur presumably over the same 
range of temperatures as the enhanced diamagnetic re- 
sponse. Both derive from the same mesocopic supercur- 
rents. The absence of significant (T2 contributions above 
Tc in the present theory is related to the fact that the 
pseudogap does not contribute to the superfluid density. 
Our derivation of the current-current correlation function 
centered on this important constraint. Thus, it should 
not be surprising that we find little signature of the su- 
perfluid density in a2 of the normal phase. Here, more- 
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over, fluctuation effects are not present since ours is ba- 
sically a mean field approach. As speculated in Ref. |31 
one should distinguish these near critical fluctuations 
from preformed pairs that persist to much higher tem- 
peratures closer to T*. 

It is notable that even cti is suppressed in the low fre- 
quency regime, when a pseudogap is present. This is 
because there are fewer fermions available to contribute 
to transport; their number is reduced since they are tied 
up into pairs. However, once the frequency is sufficiently 
high to break the pairs into individual fermions, the con- 
ductivity rises, leading [22[ to a second peak at higher 
Lj K, 2 A) in ai{uj). One can think of this two compo- 
nent behavior (which is observed experimentally [s^) as 
reflecting a transfer of spectral weight as implied by the 
conductivity f-sum rule. The behavior of CT2(i^), shown 
in Fig. mb) is rather similarly constrained by these same 
effects. On general principles, a2 must vanish at strictly 
zero frequency as long as the system is normal. Here, too, 
the low frequency behavior is suppressed by the presence 
of a pseudogap as a result of the gap-induced decrease in 
the number of carriers. The second peak in ai which re- 
flects the breaking of pairs, leads, via a Kramers Kronig 
transform to a slight depression in CT2(w) in this frequency 
regime. Thus, fT2(w) is significally reduced relative to a 
traditional Fermi liquid. 

FigElJc) compares the effects of a d-wave pseudogap on 
the diamagnetism and dc conductivity. Here we presume 
that the chemical potential of the underdoped cuprates 
is somewhat away from the Van Hove points, to avoid 
the large normal state paramagnetism [40|. The left and 
right hand axes plot the zero frequency conductivity as 
a function of varying pseudogap energy scale Apg and 
the orbital susceptibility with varying Apg respectively. 
Once pre-formed pairs are present, the diamagnetic con- 
tribution is significantly enhanced, relative to the very 



weak band diamagnetism found at IS.^ 



0. 



These observations suggest an important anti- 
correlation between the dc conductivity and the orbital 
susceptibility in the pseudogap phase, which is shown in 
FiglHfc). The dc conductivity is reduced with increas- 
ing Apg due to the opening of a gap in the fermionic 
spectrum. By contrast the existence of bosonic degrees 
of freedom (in the vicinity of condensation) and d-wave 
pairing allows for an enhanced diamagnetic response. 



CONCLUSION 

Because of the rather widespread adoption of a 

broadened BCS form for the pseudogap self energy (as 
given in Eq. ([l} ) , it seems appropriate to explore the con- 
sequences for transport, and in particular diamagnetism. 
We have done so here in a fashion which is rather inde- 
pendent of the microscopic details. The logic we followed 
is straightforward. We used the absence of a normal state 



Meissner effect to constrain the g = 0, = component 
of the current-current correlation function P (0.0). We 
extended this correlation function to finite 4- vector argu- 
ments P [Q) in a fashion which analytically satisfies the 
f-sum rules, and yields known results for limiting cases. 
This correlation function expanded for small wavevector, 
then, directly yields the orbital susceptibility, associated 
with the pseudogap self energy. 

In this way, our numerical studies are expected to 
yield rather generic results, a main summary of which 
is presented in Figure 1. The orbital susceptibility asso- 
ciated with a tightbinding bandstructure is very different 
from more familiar Landau diamagnetism [20| and, im- 
portantly, near Van Hove singularities this susceptibility 
is paramagnetic. However, the effect of a pseudogap can 
lead to large negative corrections to the bandstructure 
predictions, for the d wave case. That such diamagnetism 
is present only for d-wave symmetry is due in part to the 
presence of nodal fermions and to the extended size of 
the d-wave pairs. 

It is of interest to contrast our approach with other 
theories in the literature. Podolsky et al. ,41] and Eckl 
and Hanke [38] have respectively applied phase fluctua- 
tion approaches to address diamagnetism and conductiv- 
ity experiments leading to several predictions: (i) that 
[38| the closely related 0-2 (i^) i^"^, at sufficiently high 
w. (ii) that I4l|, the underdoped cuprates behave as a 
dilute vortex liquid over a wide range of temperatures 
above Tc. More recently, however, there has been some 
concern raised (sij about this vortex plasma model for 
the cuprates. A comparison of the conductance and or- 
bital susceptibility suggests that these vortices must ex- 
hibit an anomalously large vortex diffusion constant. We 
note that these concerns do not apply to the present ap- 
proach to transport. Also notable is a large body of work 
on related vortex scenarios 42. |43| . as well as theoreti- 
cal studies [|44] which address the superconducting fluc- 
tuation contribution to conductivity and diamagnetism 
in the presence of current-current interactions in the t-J 
model. 

The calculations we present here are different primarily 
because we consider non-Fermi liquid aspects to be dom- 
inant. We have shown how a large diamagnetic (linear) 
response derives in part from the well known enhanced 
transport contribution associated with bosons near con- 
densation. At the same time, in the present theory, these 
bosons are meta-stable at temperatures much higher than 
in the critical regime, away from where standard fluc- 
tation theories apply. This high temperature stability 
is due to stronger-than-BCS attractive interactions. We 
caution that this lattice diamagnetism is associated with 
two constraints: that the pairing be d-wave and that the 
(near mid-band) chemical potential /x lie away from the 
Van Hove singularities. When fj. is at the Van Hove point, 
the strong paramagnetism of the normal state [40j in- 
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hibits a diamagnetic response, as does s-wave pairing in in the literature so our results should have a wider appli- 
a tight binding band. cability and generality. 

In summary, this preformed pair pseudogap scenario 
leads to very different behavior for the conductivity and 

the diamagnetism and this appearsto bear on recent. This work is supported by NSF-MRSEC Grant 

otherwise challenging experiments |3ll Isg. In the for- 0820054. We thank Hao Guo and Chih-Chun Chien for 
mer case because there are fewer fermions around to valuable insights, and Dr. Guo for the analytic proof of 
contribute to the a; = conductivity (they are tied up ^j^g longitudinal f-sum rule, 
into "bosons" ) this leads to a reduction in the dissipative 
contribution to the conductivity Because we are working 
at effectively zero magnetic field, we have not addressed 
diamagnetism associated with a non-linear response, al- 
though this appears to be very anomalous experimen- 
tally 45]. We end by reiterating that our starting point. Appendix A: Effects of finite 7 
Eq. ([T|), for deriving (j{Ld) and is frequently adopted 

I 



The current-current correlation function will be discussed below in the normal state in the presence of lifetime 7 ^ , 
which is kept arbitrary in the analysis. 

The starting point in the derivation of the current-current correlation function P (q, iflm) is the phenomenological 
cuprate self-energy "^pg^K as determined by ARPES experiments. It takes the form 

^pg.K = -i7Sgnw„ -I- . (24) 

«w„ + + «7Sgnw„ 

where K = (k, ia;„) is a 4-vector, ia;„ is a fermionic Matsubara frequency, is the normal state fermion dispersion, 
7~^ is the pseudogap pair lifetime, sgn(..) denotes the sign function, and Apg ^ = '-^pg'P'k where Apg is the magnitude 
of the pseudogap and (pi^ is the gap form factor. Eql21]is the generalization of the phenomenological ^pg^K to finite 
temperature, rather than the ia;„ — >■ w -I- iO^ limit that it is normally displayed in. The other physical constraint that 
will be used is the absence of the Meissner effect above Tc- This condition is expressed 

P0,0+— =0 25 
m 

and relates the zero momentum and frequency current-current correlation function to ^/m. Eq |25l allows one to 
avoid issues of Ward identities and renormalization of the bare electromagnetic vertex when deriving , 0) because 
the tensor ^ /m is defined solely in terms of the normal state dispersion and the Green's function Gk- By 
incorporating Eg 1241 Eal25l can be used to calculate the electromagnetic vertex renormalization by first rewriting 

y y 

n jm 

^ - ^^akak^^ - ak ^ ^^^^ 

where the second equality follows from an integration by parts. The derivative of the Green's function can be written 
in terms of a bare Green's function and EalMl 




Inserting EqI27]into EqESl yields 



where Gq ^ is defined 



9k " 9k 
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and takes the appearance of a normal state Green's function with a self-energy contribution due to a scattering process 
resulting in the 17 factor. Eq l28l implies that the zero momentum and zero frequency current-current correlation 
function is 

V(0,0) = 2Y^Gj,^^ (^1 - A^kG^^ j (30) 

The key observation to make of EqllD] is that (i) the phenomenological self-energy and (ii) ensuring the absence of 
a Meissner effect above Tc were sufficient assumptions to calculate the renormalization of the bare electromagnetic 
vertex Ak in the q, i^lm = limit, namely that the dressed electromagnetic vertex Ak is 

A. = ^^Ak = ^(l-A^„kG2,v) (31) 

The transport properties derived from Eq|3I] and its extension to finite momentum q are of very general character 
in describing the pseudogap state. The remaining step is the extension of Eq l30l to finite q and iflm which will be 
accomplished at this stage by analogy to the BCS current-current correlation function but will be supplemented by 
an analytic proof of the tranverse f-sum rule that depends on the iflm dependence of P (q, ifim). One can utilize 
the q and iflm dependences in BCS theory to build further confidence by noting that the correction to the bare 
electromagnetic vertex is of the BCS-form but with opposite sign. Mirroring the q and iflm dependence of Ak, the 
Green's function, and the correction to Ak, the extension of Eg 1501 is 

P (q, inm) = 2Y^GkGk+q ^^q-^'^ ~ ^P9,k Apg,k+qG2; _kGo -k-q^ (32) 

The analogy between Eg 1321 and its BCS counterpart can be further exploited by introducing the function Fpg^K, 
defined as 

Fpg,K = ^pgM^O.-K^K (33) 

While this is of similar form to the anomalous propagator F^cK, Fpg^K does not reflect phase coherent pairs. Inserting 
Eql33]into Eql32] leads to a compact form 



K \ 



^pg,KFpg^K+Q (34) 



The proof of the transverse sum rule will be facilitated by working with the spectral representations of Gk and Fpg x, 
which are implicitly defined through the relations 

Gk = [^^£^ (35) 

/ 27r lUJn - w 



Inserting Eql35] into Eql32] leads to the result 



" ^ 9k ak j [2T,Y\uj-uj' + in„,) 

X f Ag{K w)^G(k + q, J) - Af^^ (k, lo)Af^^ (k + q, lu') 



(37) 



The current-current correlation function can be further simplified by replacing one of the spectral functions in each 
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term with the full G or Fpg function, resulting in the final expression 

PidM = 2J2^ ai^i2^-^H (38) 

k 

X ^^G(k)G(k + q, w + irj^) + AG(k + q, w)G(k, w - mm) 

-AFpg(k)i^pg(k + q,cj + - AFpg(k + q, w)i^pg(k, w - il^,„)j 

Next we prove the transverse f-sum rule in the pseudogap state for this more general case of arbitrary dissipation. 
Recall that the f-sum rule is expressed 

^^},f-(-h^'^i^,^^)]=- (39) 

where n+ = n + iO+. 
From Eq. ([28]) we have 

^ - -Eft ('^i- w 

AaiK c^)AG{k, u') - Ap^^ (k, L.)Af^^ (k, u') 
a; — w 



From Eq. we have 

k 

X ^AG(k, w) AG(k, + f]) - Ai.^^ (k, w)Af^^ (k, + f]) j 

Thus 

X {^Ag (k, a;)^G(k, c.') ~ Ap^^ (k, w) A^^,^ (k, c.') j 



(42) 



which leads to the desired result. 
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